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Abstract 

Massive rare particles have been searched for in the penetrating 
cosmic radiation using the MACRO apparatus at the Gran Sasso 
National Laboratories. Liquid scintillators, streamer tubes and nuclear 
track detectors have been used to search for magnetic monopoles 
(MMs). Based on no observation of such signals, stringent flux limits 
are established for MMs as slow as a few lO^^c. 

The methods based on the scintillator and on the nuclear track 
subdetectors were also applied to search for nuclearites. 

Preliminary results of the searches for charged Q-balls are also 
presented. 



1 INTRODUCTION 

One of the primary aims of the MACRO experiment at the Gran Sasso 
underground Laboratories is the search for magnetic monopoles at the mass 
scale of Grand Unified Theories (GUTs) of the electroweak and strong 
interactions O with a sensitivity well below the Parker bound (10^^^ cm^^ 
s~^ sr"-*^) [^ in the velocity range 4 • 10~^ < /? < 1, /3 = vjc. 

MACRO has three subdetectors: liquid scintillation counters, limited 
streamer tubes and nuclear track detectors (CR39 and Lexan) arranged in 
a modular structure of six "supermodules" (SM's). Each SM is divided into 
a lower and an upper ("Attico") part and comes with separate mechanical 
structure and electronics readout. The full detector has global dimensions 
of 76.5 X 12 X 9.3 m'^ [^ and provides a total acceptance to an isotropic 
flux of particles of ~ 10, 000 m^sr. The detector has been built and 
equipped with electronics during the years 1988-1995. Data taking began 
in 1989 with the first SM; since the fall of 1995 it is running in its final 
configuration. The response to slow and fast particles of the scintillators, 
streamer tubes and nuclear track detectors was experimentally studied 
[Q, 1^, ^ . The three subdetectors ensure redundancy of information, cross- 
checks and independent signatures for possible MM candidates. 

The analyses presented here, based on the various subdetectors in a 
stand-alone and in a combined way, refer to direct detection of bare MMs 
of one unit Dirac charge {g^ = 137/2e), catalysis cross section acat < 1 mb 
and isotropic flux (we consider MMs with enough kinetic energy to traverse 
the Earth); this last condition sets a P dependent mass threshold (~ 10^^ 
GeV for /? ~ 5 • 10~^, and lower for faster MMs). Since no MM candidate 
was found we quote new flux upper limits at the level of 2.5 • 10^^^ cm^^ 
g-i gj.-i for /? > 5 • 10^^. 

"Strange Quark Matter" (SQM) should consist of aggregates of 
comparable amounts of u, d and s quarks; it might be the ground state 
of QCD [0]. If bags of SQM were produced in a first-order phase transition 
in the early universe, they could be candidates for the Dark Matter (DM), 
and might be found in the cosmic radiation reaching the Earth. SQM in the 
cosmic radiation is commonly known as "nuclearite" and "strangelet" [^ . 

Q-balls should be aggregates of squarks, sleptons and Higgs fields |^, p^ . 
They could have been produced in the early Universe, and may contribute 
to the Cold Dark Matter. Heavy Q-balls may have originated in the course 
of a phase transition or they could have been produced via fusion processes, 
reminiscent of the big bang nucleosynthesis. Small Q-balls could also be 



pair-produced in very high energy colhsions. Rehc Q-balls can be separated 
in two classes: Supersymmetric Electrically Charged Solitons (SECS) and 
Supersymmetric Electrically Neutral Solitons (SENS). 

Some of the methods used for the MM searches may also be applied 
to search for nuclearites and for charged Q-balls (SECS). We quote upper 
limits for /3 > 5 • 10~^. 

2 SEARCHES FOR MAGNETIC MONOPOLES 

A flux of cosmic GUT supermassive magnetic monopoles may reach the 
Earth. The velocity spectrum of these MMs could be in the range 4 • 10~^ < 
P < 0.1. Our searches for MMs exploit their energy loss mechanisms in each 
of the three MACRO subdetectors. 

In scintillators the fraction of energy loss which is effective for the 
detection is the excitation energy loss which leads to the emission of light; 
in streamer tubes it is the ionization energy loss in the gas; in nuclear track 
detectors it is the Restricted Energy Loss (REL), i. e., the energy deposited 
within ~ 10 nm from the MM trajectory. In Ref. |11] a thorough analysis 



is made of these losses and of their dependence on the MM velocity. 

Independent and combined monopole analyses were performed using the 
scintillator, streamer tube and nuclear track subdetectors in different ranges 
of velocity. As already stated the results presented here apply to bare qd 
MMs and dcat < 1 mb; for the streamer tube analysis the dependence of the 
results on acat is discussed below. 

2.1 Searches with scintillators 

The searches with the liquid scintillator subdetector use different specialized 
triggers covering specific velocity regions; the searches are grouped into 
searches for low velocity (10"'* < (3 < 10~^), medium velocity (10^^ < 
13 < 10^^) and high velocity (/? > 0.1) particles. 

2.1.1 Low velocity monopole searches 

Previous searches using data collected with the Slow Monopole Trigger 
(SMT) and Waveform Digitizer (WED) were reported in Ref. ||l2|, see 
curves "A", "B" in Fig. 1. A new custom made 200 MHz WED system 
was implemented in 1995 which improves by at least a factor of two the 
sensitivity to very slow monopoles (/? ~ 10~^) and by over a factor of five 



the sensitivity to relativistic monopoles with respect to previous conditions. 
The sensitivity of the SMT/WFD was tested with LED pulses, of ~ 6.3 //s 
duration, corresponding to /3 ~ 10~^, down to the level of few tens of single 
photoelectrons, which is the signature of a slow monopole. A waveform 
analysis procedure consisted in scanning off-line the corresponding wave 
forms and in software simulation of the function of both the analog and 
digital part of the SMT circuitry on an event-by-event basis. We plan to 
report on this search in the near future. 

2.1.2 Medium and high velocity monopole searches 

The data collected by the PHRASE (Pulse Height Recorder and 
Synchronous Encoder) trigger are used to search for MMs in the range 
1.2 • 10"^ < P < 10"^ [||, |l3|. The events are selected requiring hits 
in a maximum of four adjacent scintillation counters, with a minimum 
energy deposition of 10 MeV in two different scintillator layers. Events with 
1.2 • 10~^ < /? < 5 • 10~^ are rejected because their pulse width is smaller 
than the expected counter crossing time; events with 5 • 10~^ < (3 < 10^^ 
are rejected because the light produced is much lower than that expected 
for a MM. The analyses refer to data collected by the MACRO lower part 
from October 1989 to the end of 1999 and also by the Attico from June 1995 
to the end of 1999. No candidate survives; the 90% C.L. flux upper limit is 
2.6 • 10"^^ cm-2 g-i sr~i (curve "D" in Fig. 1). 

A previous search for MMs with (3 > 10~^ based on the ERP (Energy 
Reconstruction Processor) trigger |l^, |^] is included in Fig. 1 ( curve "C"). 

2.2 Search using the streamer tubes 



The streamer tube search was described in Refs. |12, 14 1. The detection 
of MMs of 10~^ < (3 < 10"'^ is based on the Drell and Penning effects in 



the gas mixture (73% He and 27% n-pentane) filling the tubes [14, 15|. The 
analysis is based on the search for single tracks in the streamer tubes and on 
the measurement of the velocity with the "time track" . Only the horizontal 
streamer planes of the lower MACRO structure are used in the trigger; the 
Attico and the vertical planes are used for event reconstruction. Data were 
collected from January 1992 to end of 1999 for a live-time of 6.6 • 10^ hours. 
The trigger and the analysis chain were checked to be velocity independent. 
The global efficiency was estimated by computing the ratio of the rate of 



single muons reconstructed by this analysis to the expected one [14|. The 



overall efficiency was 74%. The detector acceptance, computed by a Monte 
Carlo simulation including geometrical and trigger requirements, is 4250 m^ 
sr . No monopole candidate was found. For 1.1 • 10^^ < /3 < 5 • 10^^ the 
flux upper limit is 3.1 • 10^^^ cm~^ s^^ sr~^ at 90% C.L. (Fig. 1, curve 
"Streamer" ) . 

2.2.1 Catalysis of nucleon decay 

Detailed Monte Carlo simulations were performed to study the effects of 
nucleon decay catalyzed by a GUT magnetic monopole on the streamer 
monopole trigger and its effects on the relative analysis. Both the physical 
process and the detector response were introduced in the code, taking into 
account the theoretical predictions on the cross section and on the decay 
channels. 

Many samples (each of 10,000 monopole events) were generated with 
constant monopole velocity (3 = 10~^, 5 • 10~^, 10~^, 5 • 10~^ and 10~^, and 
with catalysis cross sections acat = 10~^^, 10~^^, 5- 10~^^ and 10~^^ cm^ (in 
the last two cases other samples with (3 = 2-10^^ and 2-10^^ were produced); 
also a sample with no catalysis was simulated as a term of reference. These 
simulations were performed according to two different theoretical models 
for the catalysis cross section, one which considers it to have the same (3 
dependence for both protons and neutrons and a second one which assumes 
it to be /3-enhanced in case of protons. 

All samples were analyzed with the same program used for the real 
data. This allowed a study of the detection efficiency as a function of (3 
and (Teat- As a consequence new upper limits can be established which take 
into account also this process. Figs. 2 and 3 show for a catalysis event with 
P = 10^^ and acat = 10"^^ cm^ the time and wire views, respectively: in 
the time view the monopole straight track and the catalysis hits are clearly 
distinguishable. Fig. 4 shows the distributions of the reconstructed /?: the 
distributions are exactly peaked on the input values, which means that the 
reconstruction code works well also in the presence of catalysis hits. Finally 
Fig. 5 shows the upper limits vs /? for different acat- for low catalysis cross 
sections {(Jcat ^ 10~^^ cm^) the difference is negligible, while for higher 
values it becomes more important. 

A new analysis is in progress searching for catalysis events in the real 
data. Moreover checks are being carried out to see how the catalysis may 
affect the combined fast monopole analysis (Sect. 2.4). 



2.3 Search using the nuclear track subdetector 

The nuclear track subdetector covers a surface of 1263 m^ and the acceptance 
for fast MMs is 7100 m^ sr. The subdetector is used as a stand-alone detector 
and in a "triggered mode" by the scintillator and streamer tube systems. A 



detailed description of the method of searching for MMs is given in Ref. |16]. 
On May 2000 we began the massive etching of the CR39 sheets using the 
Bologna and Gran Sasso facilities, at the rate of about 40 m^/month. An 
area of 368 m^ of CR39 has been analysed, with an average exposure time 
of 8.5 years. No candidate was found; the 90% C.L. upper limits on the MM 
flux are at the level of 3.7 • 10"^^ cm'^ s"^ sr"^ at /3 ~ 1, and 5.4 • 10"^^ 
cm-2 s"^ sr-i at /3 ~ lO"'' (Fig. 1, curves "CR39"). 

2.4 Combined searches for fast monopoles 

A search for fast MMs with scintillator or streamer tubes is affected by the 
background due to energetic muons with large energy losses (the nuclear 
track detector is not affected). Two analyses, which combine the use of the 
three subdetector systems, were performed in order to achieve the highest 
rejection imposing looser requirements. 

2.4.1 Streamer tubes+ERP 

The analysis procedure is based on the scintillator and streamer tube data; 
the nuclear track detector is used as a final tool for rejection/confirmation of 
the selected candidates. The trigger requires at least one fired scintillation 
counter and 7 hits in the horizontal streamer planes. Candidates are selected 
on the basis of the scintillator light yield and of the digital (tracking) 
and analog (pulse charge) information from the streamer tubes. A further 
selection is then applied on the streamer tube pulse charge. After corrections 
for gain variations, geometrical and electronic non-linear effects [17|, a 90% 



efficiency cut is applied on the average streamer charge. Possible candidates 
(~ 2/year) are analysed in the corresponding nuclear track detector modules. 
The analysis refers to about 36,980 live hours with an average efficiency 
of 77%. The geometrical acceptance, computed by Monte Carlo methods, 
including the analysis requirements, is 3565 m^ sr. No candidate survives; 
the 90% C.L. flux upper limit is 6.3 • W~^'^ cm'^ g-^ sr^^ for MMs with 
5 • 10-3 < p < 0.99 (curve "E" in Fig. 1). 



2.4.2 PHRASE+Streamer tubes 

MMs with P > 10~^ are searched for by combining the streamer tube and 
PHRASE triggers. Streamer tubes are used to reconstruct the trajectory and 
pathlength, scintillators are used to measure the velocity and the light yield. 
Selected events (~ 50/year) have a single track and an energy deposition 
> 200 MeV in three scintillator layers. The event energy loss is compared to 
that expected for a monopole with the same velocity. The analysis refers to 
about 8528 live hours from May, 1997 to June, 1998. No candidate survives. 
The geometrical acceptance, including analysis cuts, is 3800 m^ sr. The 90% 
C.L flux upper limit is 2.3 • 10~^^ cm~^ s~^ sr~^ (curve "F" in Fig. 1). 

3 SEARCHES FOR NUCLEARITES 

The main energy loss mechanism for nuclearites passing through matter is 
elastic or quasi-elastic collisions Q: 

dE 2 

-— =apv 
ax 

where a is the nuclearite cross section, v its velocity and p the mass density 
of the traversed medium. 

For nuclearites with masses M > 8.4 • 10^^ GeV (~ 1.5 ng) the cross 
section may be approximated as: 



(T ~ vr 
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where pN (the density of SQM) is estimated to be pM ~ 3.5 • 10 g/cm' 
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[18|. For lighter nuclearites the collisions are governed by their electronic 
clouds, yielding o" ~ vr • 10^^^ cm^. 

Assuming galactic velocities, /3 ~ 2 • 10~^, nuclearites with masses 
< 5 • 10^^ GeV cannot reach the detector; for 5 • 10^^ < M < 10^^ GeV 
only downward going nuclearites can reach it; for M > 10^^ GeV nuclearites 
can reach MACRO from all directions. Scintillators are sensitive to the 
blackbody radiation emitted along the heated nuclearite paths down to 
/? ~ 5 • 10~^. The CR39 is sensitive to nuclearites down to /3 ~ 10~^ |19]. 
The density of the gas mixture in the streamer tubes is too low to produce 
energy losses yielding ionization, so the streamer tubes are not useful for 
nuclearite searches. 



Individual flux limits for nuclearites from the scintillator and CR39 
subdetectors are presented in Fig. 6; curves "a - d" refer to earlier searches 
with scintillators [O]; curves "e" and "f" are the updated limits obtained 
using the PHRASE system (Sect. 2.1.2) and CR39 nuclear track detectors, 
respectively. 

4 SEARCHES FOR CHARGED Q-BALLS 



The detection of Q-balls in MACRO is discussed in Ref. [20|. In that work 
it is assumed that the main contribution to the energy losses of SECS 
passing through matter with velocities in the range 10~^ < (3 < 10~^ 
are due to the interaction of the SECS positive charge with the nuclei 
(nuclear contribution) and with the electrons (electronic contribution) of 
the traversed medium. Other energy loss mechanisms, like nuclearite energy 
loss, could be considered and are under investigation. 

The MACRO subdetectors are sensitive to SECS for any value of 
the electric charge Zq > 1. The observational signatures of SECS are 
characterized by substantial energy release along a straight track with no 



attenuation throughout the detector |1C]. The energy losses of SECS in 
different subdetectors are given in Ref. [^. CR39 and scintillators are 
sensitive to Q-balls with Zg > 1 for 3 • 10^^ < /? < 0.1 and /3 > 6 • 10-^ 
respectively; in streamer tubes the detection threshold is at /3 = 2-10"'^. We 
checked that the methods to search for MMs based on the PHRASE system 
(Sect. 2.1.2 ), on the streamer tube system (Sect. 2.2) and on the nuclear 
track detector (Sect. 2.3) can be applied to SECS. Flux limits for Q-balls 
with Zq>1 obtained with the MACRO scintillators (curves "PHRASE" ) 
and with streamer tubes (curve "Streamer") are shown in Fig. 7; the limit 
obtained with CR39 (curve "CR39" in Fig. 7) applies to Q-balls with Zq = 1. 

5 CONCLUSIONS 

No MM, no nuclearite, no Q-ball candidates were found in any of these 
searches. 

The 90% C.L. flux limits for MMs versus (3 are shown in Fig. 1. The 
global MACRO limit is computed as 2.3/Xtotai where Xtotai = J2i^iJ ^^^ 
the Xl are the independent time integrated acceptances of different analyses. 
This limit is compared in Fig. 8 with the limits of other experiments which 
searched for bare MMs with g = go and (Teat < 1 mb [21-28]. 



Following the same procedure used for MMs, we obtain the 90% C.L. 
global MACRO limit for an isotropic flux of nuclearites (with masses 
> 6-1022 GeV/c2, Fig. 6); at /? = 2-10-3 the limit is 2.1-10-1^ cm^^ g-i sr^i. 
The MACRO limit for a flux of downgoing nuclearites is compared in Fig. 9 
with the limits of other experiments [25,29-31]. The Galactic Dark Matter 
(DM) limit in Fig. 9 was estimated assuming that ^max = PDMv/{27rM), 
where pdm — lO^^^i g/cm^ is the local DM density, and M and v are the 
mass and the velocity of nuclearites, respectively. 

There is a close connection between searches for Q-balls (type SECS), 
nuclearites and magnetic monopoles. The liquid scintillators, the limited 
streamer tubes and the nuclear track detector CR39 are sensitive to SECS. 
The limits for MMs obtained with the scintillators using the PHRASE 
system and with streamer tubes can be applied to SECS. The global 
MACRO limit for SECS with electric charge Zq = 1 is obtained following 
the same procedure as for MMs and nuclearites and is shown in Fig. 7 (curve 
"MACRO"). 
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Figure 1: The 90% C.L. upper limits for an isotropic flux of supermassive magnetic 
monopoles obtained using the three MACRO subdetectors: liquid scintillators 
(curves A-D), streamer tubes (curve " Streamer"), nuclear track detectors (curves 
"CR39" ) ; curves E and F were obtained combining the three subdetectors together 
(see text). The bold line is the present MACRO global limit. 
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Figure 2: Space view (XZ and YZ in the upper and lower part, respectively) of 
a simulated catalysis event. The monopolc track and the catalysis hits are clearly 
distinguishable. 
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RUN 18514 EVENT 16 B= -12 
ToF = -34241. ns 



Figure 3: Time view of the simulated event shown in Fig. 2. In the lower part 
the content of the whole 680/zs QTP memory is plotted; the upper part shows a 
magnification of the region around the fired /3 slice. The MM time track is a straight 
line; the catalysis hits are grouped in a narrow time window. 
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Figure 4: Distributions of the reconstructed /3 for 4 simulated catalysis cross 
sections. The distributions are exactly peaked around the input values, which 
means that the recontruction procedure (used in the standard streamer analysis), 
gives the correct (3 even in presence of catalysis hits. 
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Figure 5: Flux limits from streamer tubes for 4 values of acat, which is assumed to 
be the same on protons and neutrons: (from above) acat = lO"^**, 5 • 10~^^, 10~^^ 
and 10~^^ cm^. The simulation was performed at fixed /3; the lines are only a guide 
to the eye. 
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Figure 6: The 90% C.L. upper liniits for an isotropic flux of nuclearites obtained 
using the liquid scintillator (curves "a" - "e") and the CR39 nuclear track (curve 
"f") subdetectors; the bold line is the MACRO global limit. 
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Figure 7: The 90% C.L. upper limits for an isotropic flux of charged Q-balls 
obtained using the liquid scintillator (curve "PHRASE"), the CR39 nuclear track 
detectors (curve "CR39") and the streamer tubes (curve " Streamer"). The bold 
line is the MACRO global hmit. 
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Figure 8: The global MACRO 90% C.L. upper limit for an isotropic flux of g=gD 
magnetic monopoles extends from /3 — 4 ■ 10^^ to 1; it is compared with the limits 
obtained by other experiments; at values of /? ~ 1 we show the limits from the 
Baikal and Amanda collaborations. 



19 



7 

CM 
I 

E 
o 

_i 

6 

K 
o 

CD 



0) 

Q. 

D 
X 




10 



24 



Nuclearlie m.ass (Ge\//c^) 



Figure 9: The global MACRO 90% C.L. upper flux limit for nuclearites with 
j3 — 2 ■ lO^'' at ground level, versus nuclearite mass, is compared with the limit 
obtained by other experiments and with the galactic DM bound. The limit above 
Mat > 6 • 10^^ GeV corresponds to an isotropic flux; for Mn < 6 • 10^^ GeV the 
limit corresponds to nuclearites reaching the detectors only from above. 
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